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Abstract 

For the fully anisotropic simple-cubic Ising lattice, the critical finite-size 

scaling amplitudes of both the spin-spin and energy-energy inverse correlation 

lengths and the singular part of the reduced free-energy density are calculated 

by the transfer-matrix method and a finite-size scaling for cyclic L x L x oo 

clusters with L = 3 and 4. Analysis of the data obtained shows that the ratios 

and the directional geometric means of above amplitudes are universal. 
PACS numbers: 05.05.+q, 75.10.H, ll.15.Ha 
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I. INTRODUCTION 



Conforming with the Privman-Fisher hyperuniversahty hypothesis, the finite-size scahng 
(FSS) equations for the inverse correlation lengths and the singular part of the reduced free- 
energy density near the bulk phase transition of a system have, respectively, the form (for 
reviews see Refs. 

^,,L{t, h) = L-^X,{CitLy^, C2hLy^) (1) 

and 

fi'\t, h) = L-^Y{C,tLy^, C2hLy^) . (2) 

Here L is a characteristic size of finite or partly finite subsystem, the index i labels the 
types of correlation lengths [spin-spin (z = 1), energy-energy {i = 2), etc], d is the space 
dimensionality, t = (T — Tc)/Tc, h is an external field, yx and i/h are the critical exponents, 
Xi{x, y) and Y{x, y) are the scaling functions which, within the limits of universality classes, 
can else depend on the type of boundary conditions and the subsystem shape; all non- 
universality of a model is absorbed in the metric factors Ci and C2. Equations (1) and (2) 
allow to find the universal combinations for the FSS amplitudes at the phase-transition point 
t = h = 0. In particular, the amplitudes for the inverse correlation lengths Ag = Xi(0,0) 
and Ae = X2(0,0) and for the free energy Af = y(0,0) must be universal themselves. In 
the case of strips with periodic boundary conditions, they are (see, e. g., Ref. 0) 

vrc 

As = 7rr], Ae = 7n]e, Af = (3) 

where t] and rj^ are the exponents of the decay law correspondingly of the spin-spin and 
energy-energy correlation functions {1] = 1/4 and rje = 2 for the fiat Ising model) and c is 
the central charge of Virasoro algebra (c = 1/2 for the two-dimensional Ising lattice). 

All foregoing statements are applied to the spatially isotropic systems. Lattice anisotropy 
is a marginal effect and hence the amplitudes and their combinations, strictly speaking, must 
depend on anisotropy parameters . However in the case of the anisotropic two-dimensional 
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Ising model, it has been established that although the inverse correlation-lengths and 
free-energy amplitudes get a non-universal factor, {a labels the directions along which 
an L X cxD strip is infinite; here, a = x,z), it is common and the directional geometric mean 
R = {RxRzY^'^ is a constant (equalling the unity). Therefore, the universality is preserved 
for the ratios and the directional geometric means of these amplitudes. 

In the light of above, it would be interesting to clear the matter up in three dimensions. 
Such attempt is undertaken in the present paper. We consider the three-dimensional Ising 
model on a simple-cubic lattice with different interaction constants Jx, Jy, and Jz along all 
three spatial directions. The lattice is approximated by the L x L x oo bars with periodic 
boundary conditions in both transverse directions. Such boundaries eliminate undesirable 
surface effects and hence improve a quality of approximation. By the transfer matrix (TM) 
method combined with FSS analysis for the subsystems with sizes L = 3 and 4, we determine 
at first the critical temperatures depending upon anisotropy parameters Jx/Jz and Jy/Jx- 
(We consider a system at least with two non-zero couplings; unless otherwise stated, the 
L X L X oo parallelepipeds are taken infinitely long in the z direction.) After this, the 
FSS amplitudes of the inverse correlation lengths and the free energy are calculated at the 
critical points found. The obtained results demonstrate the independence of the amplitude 
ratios on the parameter Jx/ Jz when Jy/ Jx is fixed. Moreover, the analysis shows that the 
ratios are also independent on the second anisotropy parameter Jy/ Jx at any rate in the 
region Jy/Jx — 1- Finally, our calculations give evidence in the constancy of the directional 
geometric mean of the spin-spin inverse correlation length amplitude in three dimensions. 
Together with an invariance of the ratios, this implies that the directional geometric means 
of other amplitudes must be universal also. 



II. MODEL AND SOLUTION OF THE EIGENPROBLEMS 

The Hamiltonian of Ising model on a simple-cubic lattice with nearest-neighbor interac- 
tions reads 
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(4) 



The spin- field variables Sijk are located in the lattice sites and take the values ±1. 
The transfer matrix ^ of an L x L x cxd subsystem is introduced by elements 

L 



where = Ja/ksT (now a = x,y,z); Sn+i = Sn and = Sij by all i, j = 1,2, . . . , L. 

The matrix V is real, symmetric and has an order of 2^ where N = L"^ equals the number 
of chains in a system; that is dense and all its elements are positive. 

The principal task is to find the eigenvalues of V because, for example, the density of a 
free energy measured in units of —ksT is given by 



where Aq is the largest eigenvalue of a TM. The inverse longitudinal correlation lengths 
(mass gaps) equal 



where Ai, A2 . . . are the next (after Aq) dominant eigenvalues of TM for the subsystem. 

In order to solve the TM eigenproblem for L as large as possible, we reduce the TMs to the 
block-diagonal forms using a symmetry under the transformations of the group Z2 x T A C2V 
Here Z2 is a group of global spin inversions S ^ —S, T is a group of translations in the 
transverse directions of a bar, and C2V is the point group consisting of rotations around the 
axis of a subsystem at angles multiple to vr and the reflections in planes going through this 
axis and the middles of opposite sides of an L x L x 00 parallelepiped. 

There is no necessity to perform the full quasidiagonalization of TMs because the leading 
eigenvalues are distributed only among two subblocks. Owing to the Perron theorem 0, 
Aq lies in the subblock of an identity irreducible representation. Ai is located in the other 




(5) 



(6) 



i^i,L = ln(Ao/Ai) , 



(7) 



4 



subblock — it is built on the basis functions which are symmetrical under all transformations 
of the space subgroup T A and antisymmetrical under the transformations including a 
spin inversion. A2 is situated again in the subblock of an identity irreducible representation. 
(In connection with this see, for example, Ref. 0.) 

As a group-theoretical analysis shows (see Appendix A), both subblocks containing the 
largest eigenvalues have sizes of 18 x 18 in the case of 3 x 3 x 00 cluster. For a cylinder 
4 X 4 X 00, the TM 65 536 by 65 536 is reduced to a block-diagonal form in which the required 
subblocks have the orders 787 and 672. The final extraction of needed eigenvalues of TMs 
was carried out by a numerical solution of eigenproblems for the corresponding subblocks. 
By this, we applied the conjugate gradient method and, if necessary, used also the library 
functions tred2 and tqli [0]. Calculations were run on IBM PC-486 computer in the operating 
system LINUX. 

III. CALCULATION OF THE CRITICAL AMPLITUDES 

So, the FSS amplitudes for the inverse correlation lengths of the spin-spin and energy- 
energy correlation functions are equal to 

As = Lki^l (8) 

and 

Ae = Lk2,l , (9) 

where ki ^ and k,2,l have been taken at the phase-transition point (by h = 0). This point 
itself was determined from the renormalization-group equation 

Lm^iiTc) = (L - 1)k,^l-i{Tc) (10) 

with L = 4. The amplitude for the singular part of a free-energy density, is found from 
a system of equations 
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/l = /o + L-'Af (11) 

with L = 3 and 4. Here /o denotes the regular (background) part of a free-energy density; 
/s and /4 are taken again at the critical points. 

The critical temperatures, amplitudes, and background /o calculated at different values 
of the anisotropy parameters Jx/ Jz and Jy/ Jx are collected in table I. In Eqs. (11), the spatial 
dimensionality has been put d = 2 ioi Jy = and d = 3 for Jy ^ 0. It should be noted also 
that, for finite L, Eq. (7) leads to the wrong values for k,2,l in the limit of non-interacting 
strips {Jy = 0). Due to > A0A2 (Aq > Ai > A2 are the largest eigenvalues of a transfer 
matrix for the strip), A2 = (AqAi)^ by L = 4 and therefore K2,l = ln[AQ/(AoAi)^] = 2ki i. 
However, the correct values are given by formula K2^l = ln(Ao/A2) which has been used to 
build up the table I. 

In table II, we present the data for the directional geometric mean of the spin-spin inverse 
correlation length amplitude As. Calculations were performed by the equation 

f (4-) 4^)) tf Jy = 

^ = ' , (12) 

[ (4-)4^)4-))V3 ^f Jy^O 

where A^f'' is the amplitude of the spin-spin inverse correlation length when the bar LxLx 00 
was stretched (for given Jx, Jy, and J^) along the a direction. 

IV. DISCUSSION 

Consider first the behavior of absolute amplitudes. For the three-dimensional systems, 
available information about them is very scanty. In the periodic cylinder geometry, it seems 
to be known only the estimates for the correlation-length amplitudes found by Monte Carlo 
simulations on the fully isotropic {Jx = Jy = Jz) lattices Lx Lx 128 with L = 4,6, 8, and 10 
(Ref. 0). For the inverse correlation- length amplitudes, these estimates {L = 10) yield Ag = 
1.342 and Ae = 4.78. Appealing to table I, one can convince oneself that our calculations 
conform with these values. Note also that the available high-temperature series for the free 
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energy of a fully isotropic simple-cubic Ising lattice yields /o = 0.77711 at criticality. Our 
estimate for the background, 0.773, is in good agreement with this magnitude. 

In the two-dimensional case {Jy = 0), there exists, vice versa, complete information 
concerning the FSS amplitudes for the inverse correlation lengths and the free energy in the 
rectangular lattice with arbitrary anisotropy P,H]: 



A.. 



TT 
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A, = 27V 



and 



Af = — 
^ 12 



sinh.{2 J^/ksTc 
sinh(2 J^/fc^Tc 



sinh(2J^.//i;Brc 
sinh.{2 Jz/ksTc 



sinh.{2 Jx/ksTc 



.1/2 



1/2 



(13) 



(14) 



nl/2 



smh{2 J J kBTc) 
where the critical temperature satisfies to the equation 



(15) 



smh - — — smh 



1 . 



(16) 



^ksTcJ \kBTc, 

Our numerical results reproduce these rigorous dependencies with acceptable accuracy. For 



(17) 



the isotropic square Ising lattice, the critical free energy is (see Ref. [|T0| 



/o = 2G/n + i In 2 = 0.929 695 



((7 = 1^^ — 3^^ + 5^^ — . . . is Catalan's constant). Appropriate value from table I (/o at 
Jx = Jz and Jy = 0) agrees to within 1.7% with the given exact quantity. 

Inspecting table I, we see the amplitudes vary in wide limits reaching several orders. The 
behavior is changed into a contrary one for their ratios. First what draws attention is that 
the ratios Af./As and Af/A^ stay practically unchanged with variation of Jx/Jz on three 
orders (1 — 10~^) by given Jy/Jx- In the two-dimensional space {Jy = 0), the mean (here 
and below, over Jx/Jz) value of Aj/Ag equalling to 0.331 conforms with the true value 1/3; 
the mean of Af,/As equals to 7.2 that agrees, in order of magnitude, with the exact value, 
8, for the Ae/Ag [see Eqs. (13) - (15)]. For the three-dimensional lattice with Jx = Jy, the 



mean value of Af/Ag is 0.288. This quantity agrees with estimate Af/Ag = 0.272 which 
follows from the calculations of relative amplitudes for the inverse correlation lengths and 
the free energy in the Hamiltonian limit of a three-dimensional Ising model (square lattices 
L X L with sizes L up to 5) [|ri|]. According to table I, the ratio for the inverse correlation- 
length amplitudes is A^/Ag = 3.53(6) in the discussed case. This estimate is in agreement 
with the mean values Ae/Ag = 3.62(7), Ref. fll, and A^/Ag = 3.7(1), Ref. ||]. Thus, the 



amplitude ratios A^/Ag and Af/Ag are not only universal with respect to the Jx/Jz but also 
their values agree quantitatively with available estimates in two limited cases: Jy/Jx = 
and 1. 

We now discuss the dependence on Jy/ Jx in the intermediate region. In the limit Jy/ Jx — >■ 
0, the L X L X (yo bar decomposes into L of independent strips L x oo and consequently the 
TM of the bar is factorized into the direct product of TMs for the strips. Since the TM of 
the bar is finite by finite L, its eigenvalues are continuous functions of model parameters. 
Hence there must exist the d = 3 —>■ d = 2 transition region when Jy/Jx ^ 0. To estimate 
its sizes by using L, we have calculated the critical exponents u and 7/1/. The calculation 
was performed via the ordinary FSS formulae (see, e. g., Ref. jr5[): 

ln[L/(L - 1)] 



and 



(18) 



Hxl/xl-i) . . 

^' MV(^-i)]' ^ ^ 

in which we put L = 4. Here k'^ ^ is the derivative of ki ^ with respect to the temperature 
and ^ind xl are the magnetic susceptibilities of subsystems at the phase-transition 

point. (Formulae for the susceptibilities are derived in Appendix B.) How the calculation 
gives the critical exponents v and, especially, 7/z/ are practically constants with respect 
to Jx/Jz (= 1 — 10^'^). Their dependences on Jy/Jx are shown in Fig. 1. Within the 
section 0.2 < Jy/Jx < 1, the exponents u and 7/1/ preserve the unchanged values equalling, 
respectively, to 0.67 and 1.97 that agrees with available estimates for these exponents in the 



case of the fully isotropic three-dimensional Ising model (Ref. and references therein). By 
Jy = 0, our calculation yields u = 1.06 and ^jv = 1.74. These magnitudes conform closely 
with the exact values of discussed exponents in two dimensions: v = \ and ^jv = 7/4. In 
Fig. 1, it is clear-cut displayed the region < Jyj Jx < 0.1 — 0.2 where a smooth transition 
occurs from the d = 3 exponent values to the d = 2 ones. Consequently, one does not 
consider the L x L x (yo lattice with L < 4 as a three-dimensional one when Jy/Jx < 0.2. 
In order to support this conclusion, we have calculated the "effective" lattice dimensionality 
solving the system of Eqs. (11) with L = 2,3, and 4 and treating d in it as an unknown 
continuous variable d*. (For the fully anisotropic 2 x 2 x oo Ising lattice, there is an exact 
analytical solution ||15|.) The conclusion is d* does not depend on Jxl-Jz and its plot on 
Jy/Jx is also presented in Fig. 1. This plot has a more qualitative character because in 
the calculation a cluster with an extremely small size L = 2 has been used. Nevertheless, 
the presented dependence indicates that the lattice dimensionality d* is less than three by 

Jy/Jx < 0.3. 

As mentioned in Sec. Ill, the energy-energy inverse correlation length k,2,l (and hence 
the amplitude Ae) has a false behavior in the limit Jy/ Jx ^ due to finite sizes L. The 
scaling amplitude Af obtained from Eqs. (11) with L = 3 and 4 suffers from a similar 
defect. By finding of Ae and Af, it is not allowed to change the order of the limits L — >■ 
oo and Jy/ Jx 0. (Note in passing that the calculation of hi^l and Ag is free upon 
such requirement.) Taking into account these circumstances, let us consider in Fig. 2 the 
obtained dependencies of ratios A^/Ag and AeAj/A^. The plots of both dependencies have 
the horizontal sections by small deviations of Jy/ Jx from unity. Thus, the amplitude ratios 
do not depend on the second anisotropy parameter Jy/Jx in this region of its values. As 
Jy/ Jx is decreased, both quantities tend to the incorrect limits. 

The recognized properties of the critical FSS amplitudes by a given orientation (a) of an 
L X L X oo bar in the anisotropic Ising lattice can be described by equations 

«i,L,a(0,0) =L-ii?„X,(0,0) (20) 
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and 

fi'lM^L-'R^Y (0,0), (21) 

where Xj(0, 0) and Y{0, 0) are amplitudes of the isotropic model and Ra — Ra{Jx/ Jz, Jy/ Jx)- 
The given equations are true at Jy/ Jx — and, according to the presented data, when 
Jyl Jx ~^ ^- Equations (20) and (21) are likely to be vahd also over the wider range of Jy/ Jx- 
This is confirmed qualitatively by the calculation of ^/ from Eqs. (11) with L = 2, 3, and 4 
(without supposition that d = 3 for all Jy ^ 0). 

Discuss now the behavior of the directional geometric mean of the spin-spin inverse 
correlation length amplitude (table II). In the two-dimensional case (column with Jy/Jx — 
0), Ag loses a stabihty when Jx/Jz < 10~^. This is obviously connected with small widths 
of strips by which we approximate the system. Situation is perceptibly better in three 
dimensions. Here Ag — 1.7(3), i. e. the percentage error equals 18%. With such accuracy, 
we may consider As as a constant. 

V. CONCLUSIONS 

In this paper, the TM-FSS calculations of critical temperatures, exponents, amplitudes, 
and free-energy background for the fully anisotropic three-dimensional Ising model have 
been carried out. The data obtained allow to make the following inference concerning the 
structure of critical FSS amplitudes of the inverse correlation lengths and the free energy: 
Similarly to the two-dimensional case, all lattice-anisotropy parameters are absorbed in a 
separate prefactor which is common for named amplitudes and the directional geometric 
mean of which is the unity. 
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APPENDIX A: QUASIDIAGONALIZATION OF THE TRANSFER MATRICES 

The group Z2 x T A C2V has an order g = 8L^. Its generating elements are a spin 
inversion /, translations on one step and ty, and reflections in the symmetry planes Uy 
and (Ty. In the transfer matrix space \Sii, S12, . . . , Sll), they are defined as 

-'^I'S'ii, 5'i2, . . . , Sll) — \ — Su, —S12, ■ ■ ■ , —Sll) , (Al) 

<S'i2, . . . , Sil', S21, S22, ■ ■ ■ , S2L', ■ ■ ■ ', Sli, Sl2, ■ ■ ■ , Sll) 

— \SiL, Sll, ■ ■ ■ , SiL-i', S2L, S21, . . . , S2L-1', Sll, Sli, ■ ■ ■ , Sll-i) (-^2) 

tylSu, Si2, ■ ■ ■ , Sil', S21, S22, ■ ■ ■ , S2L', ■ ■ ■', Sli, Sl2, ■ ■ ■ , Sll) 

— I'S'21, >S'22, . . . , S2L', • • • ; 'S'li, Sl2, ■ ■ ■ , Sll', Su, S12, ■ ■ , Sil) (^3) 

(^vlSu, Si2, ■ ■ ■ , Sil', S21, S22, ■ ■ ■ , S2L', ■ ■ ■ ', Sli, Sl2, ■ ■ ■ , Sll) 

— \Sli, Sl2, ■ ■ ■ , Sll', ■ ■ ■ ', S21, S22, ■ ■ ■ , S2L', Su, S12, ■ ■ ■ , Sil) (^4) 

(^vlSn, S12, ■ ■ ■ , Sil', S21, S22, ■ ■ ■ , S2L', ■ ■ ■ ', Sli, Sl2, ■ ■ ■ , Sll) 

= \SiL, ■ ■ ■ , S12, Sii; S2L, ■ ■ ■ , S22, S21', Sll, ■ ■ ■ , Sl2, Sli) ■ (^5) 

Other transformations of the group are the corresponding combinations of above operations. 
Multiplying from the left the equations like (Al)-(A5) on conjugate vectors and taking into 
account the orthonormality condition 

(-S"!!, -S'12, . . . , Sll\S[i, S[2, S'll) = (^Sii5i/5i25;2 ■ ■ ■ ^SllS'^^^^ (A6) 
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i^ss' = ||5' + 5"| is a Kronecker symbol), we find the original representation T of the group. 

All matrices of representation built commute with V. For instance, using Eqs. (5) and 
(Al), we have 

{Sii, 5*12, . . . , Sll\I -^VIlS'i-^^, 5(2, • • • , S'^^) 

= (—'S'li, —512, • • • 5 — 'S'llIV^I — S[^, —S'i2, • • • , ^S'j^l) 

= (5*11, 5*12, . . . , SLL\y\S[i, S[2, • • • , S^l) (A7) 

so that [V, /] = 0. The same is valid for all other transformations of the group. 

The traces of matrices built are characters of representation F. For the 3 x 3 x cxd case, the 
characters of original representation together with characters of irreducible representations 
F*^^^ and F*^^-* to which correspond the subblocks containing the largest eigenvalues are given 
in table III. Using this table and utilizing the formula for counting the multiplicities with 
which a given irreducible representation enters into an original representation (see, e. g.. 



Ref. [|T6l) 



a^ = -T.9^^'^*X^ (A8) 

9 i 

{gi is a number of elements in ith class, xl^"* is a character of element from ith class in 
/ith irreducible representation, and Xi is a character of element from ith class in an original 
representation) we find the composition of representation F: 

F = 18(F(i) + F(2)) + . . . . (A9) 

It follows from here that in a basis where the representation F is completely reducible the 
transfer matrix of 512-th order will take a quasidiagonal form in which both subblocks 
corresponding to the one-dimensional irreducible representations F^^^ and F*^^) will have the 
sizes 18 by 18. 

The basis vectors of irreducible representations on which the transfer matrix takes the 
discussed block-diagonal form are built with a help of projection operators In the case 
of an L = 3 subsystem, the basis vectors for the irreducible representations F^^'^^ are 
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(1.2) 


= (mi ± M5i2)/V2 


(1.2) 


= E'i ^^(ms ±M505)/V6 


(1.2) 


= E'i Gi{u74 ± M439)/V6 


(1.2) 
¥^4 


= E'.G,K5±M428)/2v^ 


(1.2) 


= E'i G,(m2±M51i)/3\/2 


(1.2) 
^6 


= E'i Gi{uA ± M509)/3V2 


(1.2) 


= E'iG'i(nio±U503)/3V2 


(1.2) 
^8 


= E'.G',(n28±U485)/3y2 


(1.2) 


= E'iGi(u79±ti434)/3 V2 


(1.2) 


= E'i ± ^i502)/6 




= Gi{ui^ ± U498)/6 


(1.2) 
Vl2 


= E'j Gi{u75 ± M438)/6 




= E'iG'i(Mi6 ±M497)/6 


(1.2) 
¥^14 


= E'i ^^(Mre ±M437)/6 


(1.2) 
¥^15 


= E'i ^^(Mao ± M483)/6 


(1.2) 
¥^16 


= E'i G'i(M84 ±M429)/6 


(1.2) 


= E'^G',(ni2±U5oi)/6v^ 


(1.2) 
¥^18 


= E'i ^^(use ± n427)/6v^ , 



(AlO) 



(All) 



where 

«i = |1, 1, 1; 1, 1, 1; 1, 1, 1), M2 = |1, 1, 1; 1, 1, 1; 1, 1, -1), . . . , 
-"512 = I - 1, -1, -1; -1, -1, -1; -1, -1, -1) • 

The plus and minus signs correspond to the basis vectors of irreducible representations T^^^ 
and r*^^\ respectively. For shortening of a listing, only the /-conjugated pairs of generating 
orths are shown in Eqs. (AlO). The numbers of orths in a pair {n and n') are connected by 
a relation ra' = 2^ + 1 — n. Acting on such orths by operators Gj G T A C2v and taking on 
each step only the new w-orths (this peculiarity is marked by prime on the sum symbol) , we 
obtain the expressions for the basis functions in explicit form. 

Finally, having the basis vectors for the irreducible representations, one can find the 
matrix elements of subblocks with the transfer matrix eigenvalues under search. For the 3 x 
3 X oo task, the matrix elements of subblocks corresponding to the irreducible representations 
r*^^'^-* have been given with all necessary coefficients in Ref. 

In the case of 4 x 4 x oo subsystem, the basis vectors of F*^^) and F*^^-* can be taken in 
the form 
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V'f ''^ = (ei ± e65536)/\/2 V'^'^ = E'^ G,(e2 ± e65535)/4\/2 

V'Sf = E'i ^,(613 670 ± e51867)/4v^ = E'^ ^,(613 674 ± e51863)/8^ (Al2) 

^673 = \ E'i Gie256 ■ ■ ■ -0787 = (623131 + 642406)7 ^2 , 

where 

ei = |1, 1, ... , 1), 62 = |1, 1, ... , -1), . . . , 665536 = | " 1", 1, ■ ■ ■ , "1) ■ (Al3) 

The basis functions (A12) from 1 to 672 and then from 673 to 787 are ordered with the 
numbers of the first generating e-orths increasing. Using Eqs. (5), (A12) and (A13), we 
evaluate the matrix elements Vil"^^ — '4>f''^^^V'4>f''^^ for subblocks corresponding to the 
irreducible representations F*^^'^) . The matrix elements are 

exp[|(m» + rrf^K^ + \{m\ + m\)Ky\ 

(A14) 

and 

exp[|(m^ + m^^K^ + \{m\ + m\)Ky\ , (A15) 

where rii are lengths of basis vectors, m" and are the reduced partial energies of spin 
configurations in orths of ith vector. All coefficients 51^*'-'^ are non-negative and satisfy to 
the "sum rules" 

gj}''^ + 2 ^ g^'^^ = minK, n,) . (A16) 

s=l 

We did not keep the coefficients g^'''^ but restored them for each matrix element V^^p from 
Eqs. (A16). As a calculation shows, the coefficients g^^'^^ with s ^ are not greater than 
60. Hence, it is enough to take one byte for every element of the g^-array, i. e. to use the 
data type 'char' in C code. Thus, it is required 2 480 624 bytes of a memory to store the 



(1) ^ max(ni,nj^ 

ij ~ 



^f^') + 2^yPcosh(2sX,) 

s=l 



(2) _ ^ max(ni,nj) 



5:^(^'^)sinh(2si^,) 

.s=l 
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(^-coefficients for a triangle part of symmetric matrix V^^\ The values of coefficients lie 
in the range from —28 to +40 and we alloted in addition the 1 809 024 bytes of a memory 
for the ^-coefficients of matrix V^'^\ 



APPENDIX B: FORMULAE FOR THE CALCULATION OF SUSCEPTIBILITIES 

In deriving of formulae for xl, we will point out from a fluctuation-dissipation relation 
connecting the susceptibility with a magnetic moment Ai (see,for example, Ref. ||T7|): 

Here ^A = Y^ijk ^ij where Sfj = Sijk is the total magnetic moment of L x L x M periodic 
subsystem; the brackets refer to average on Gibbs distribution. Taking into account the 
translational invariance of a cluster in the longitudinal (z) direction, one can write Eq. (Bl) 
in the form 

■1 M-l 

xl{t) = —— hm y: ((^11 + ^1^2 + • • • + Si,) {s',r + sir +■■■+ stti) ■ m 

To calculate the statistical means, we use the transfer matrix technique. Let us introduce 
in addition the spin matrices making by this the one-dimensional order of pair of indexes 
hj ^ 1 = L{i - 1) + j: 

5"; = 1 X . . . X 1 xcr^ X 1 X . . . X 1 , (B3) 

" T^l ' N-l 

where 1 denotes the unit matrix of second order and cr^ is Pauli's 2;-matrix; = L^. This 

allows to rewrite Eq. (B2) as 

1 1 ^^-1 . . . 

^'^^^ = itL TYV^ + • • • + + • • • + Sn)V''-^] . (B4) 

From here, by passing under trace symbol into diagonal representation of the transfer matrix 
and by taking into account the non- degeneracy of its largest eigenvalue, we obtain 

1 2^-1 A I A 
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where Fq, Fi, . . . are eigenvectors of matrix V corresponding to its eigenvalues Aq, Ai, 

Further, the operator S — Si + . . . + Spf is invariant with respect to all purely spatial 
transformations and breaks the Z2 symmetry. Therefore, the matrix elements entering into 
Eq. (B5) are not zero only for "transitions" from the identity irreducible representation F^^^ 
just into the irreducible representation F^^^. 

Vector Fq is a linear combination of basis functions only of the identity irreducible rep- 
resentation. In the case of 3 x 3 x 00, 
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p _ V- .(0) (1) 



(B6) 



i=l 



where //°'' are components of eigenvector answering to the largest eigenvalue (Aq) of subblock 
of the identity irreducible representation. Using Eqs. (AlO), we find that 



(B7) 



where 



rrii 



{9,3,3,3,7,5,5,1,1,5,3,3,1,1,1,1,3,1}; 



(B8) 



rrii is the magnetic moment of spin configurations in the ith basis vector. As a result, we 
obtain from Eq. (B5) the following work formula for a calculation of the susceptibility: 



1 18 A _L a(2) 



QksT Ao - A 



(2) 
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(0) f{i) 

j 



(B9) 



Here are components of ith eigenvector corresponding to eigenvalue A^^^* for the subblock 
of irreducible representation F*^^) . 

Analogous formula take place for the L — A subsystem: 



X4(T) = 



gAn + A 



(2) 



672 

E "..-/)"'/ 



(0) f{i) 

j 



(BIO) 



All quantities entering in this expression should be taken, of course, for the 4 x 4 x 00 model. 

Therefore, the calculation of susceptibilities requires the solution of a part eigenproblem 
for the subblock of an identity irreducible representation and the solution of a full eigen- 
problem for a second subblock which corresponds to the irreducible representation F*^^) . The 
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part eigenproblem was solved again by the conjugate gradient method and the full 
by using the library C pair tred2 - tqli. 
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TABLES 

TABLE I. Critical temperatures, background, critical FSS amplitudes and their ratios for 

different values of the anisotropy parameters Jx/Jz and Jy/Jx- 





Jx/Jz 


ksTc/Jz 


As 


Ae 


Ae/Ag 


Af 


/o 


Af/As 


1.0 


1.0 


4.58104 


1.4401 


4.9627 


3.44 


0.4189 


0.773 


0.290 




0.1 


1.35037 


0.3613 


1.2847 


3.55 


0.1044 


0.959 


0.288 




0.01 


0.65458 


0.0723 


0.2576 


3.56 


0.0208 


1.576 


0.287 




U.UUi 


n /inoi 7 


n ni 1 c; 

U.UiiO 


0/111 


O.O t 


n nn^*? 

U.UUOO 




U.ZoO 


0.75 


1.0 


4.18009 


1.3345 


4.5795 


3.43 


0.3953 


0.775 


0.296 




0.1 


1.27931 


0.3317 


1.1712 


3.53 


0.0973 


0.985 


0.293 




0.01 


0.63312 


0.0653 


0.2310 


3.53 


0.0191 


1.623 


0.292 




n nni 


u.oyyoo 


n ni nQ 

U.UiUo 


n nQf?c; 

U.UoDO 


0.04 


n nn^n 

U.UUoU 


/.OUo 


n 001 
u.zyi 


0.5 


1.0 


3.73973 


1.2288 


4.0812 


3.32 


0.3924 


0.782 


0.319 




0.1 


1.19903 


0.3005 


1.0251 


3.41 


0.0943 


1.019 


0.313 




0.01 


0.60815 


0.0580 


0.1981 


3.41 


0.0181 


1.683 


0.312 




n nni 

U.UUi 


U.ooooZ 


n nnon 
u.uuyu 


n n^no 
u.uouy 




n nn98 

U.UUZo 


9 1^78 
Z.O / o 


n ^1 1 


0.25 


1.0 


3.22427 


1.1256 


3.3542 


2.97 


0.4407 


0.803 


0.391 




0.1 


1.10117 


0.2665 


0.8151 


3.05 


0.1016 


1.073 


0.381 




0.01 


0.57655 


0.0500 


0.1533 


3.06 


0.0190 


1.767 


0.380 




0.001 


0.37453 


0.0076 


0.0235 


3.09 


0.0029 


2.675 


0.381 


0.0 


1.0 


2.32081 


0.8917 


5.9901 


6.71 


0.2952 


0.914 


0.331 




0.1 


0.91079 


0.1856 


1.3661 


7.36 


0.0616 


1.232 


0.331 




0.01 


0.51058 


0.0327 


0.2418 


7.39 


0.0108 


1.983 


0.330 
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0.001 0.34346 0.0048 0.0349 7.27 0.0016 2.915 0.333 
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TABLE II. Directional geometric mean of the spin-spin inverse correlation length amplitude 
Ag by different values of Jx/ Jz and Jy/ Jx- 



0.25 0.5 0.75 1.0 

1.0 0.891 1.57 1.46 1.43 1.44 

0.1 0.833 2.00 1.87 1.80 1.76 

0.01 0.577 2.02 2.01 2.01 2.02 

0.001 0.294 1.48 1.51 1.54 1.57 
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TABLE III. Characters of the group Z2 x T A C2v in the case L = 3; here T A C2v ^ x C3V . 



Atxty ^^^x^y-} ^^^x^vt ^^^y^v 

111 11 1 

r(2) 111 11-1 

r 512 64 32 8 4 



23 



FIGURES 

FIG. 1. Critical exponents v and ^jv (left scale) and the effective lattice dimensionality d* 
(right scale) versus anisotropy parameter Jy/Jx- 

FIG. 2. The amplitude ratios A(,/As and A^Af/A^ against the anisotropy parameter Jy/Jx- 
The curve parts which are considered as non-physical ones are shown by dashed line. 
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